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Method of Tracking Dynamic Communication Parameters during a Link Interruption 



(57) In a mobile communication system (90) wherein communication parameters vary appreciably between 
transmissions, subscriber units (200) are susceptible to link blockages resulting in a loss of successive 
feedback instructing subscriber units (200) of adjustments to communication parameters to be employed in a 
subsequent transmission on a communication link (105). When a mobile subscriber unit (200) is airborne, an 
outage may occur from a control surface blockage during maneuvering. Aircraft avionics (290) supplies 
acceleration data representing aircraft movement for use in predicting subsequent communication 
parameters. These predicted parameters are employed in re-establishing communications upon the removal 
of such a link blockage. 
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METHOD AND SYSTEM OF TRACKING DYNAMIC COMMUNICATION 
PARAMETERS DURING A LINK INTERRUPTION 

RELATED INVENTIONS 

The present invention is related to the following 
invention which is assigned to the same assignee as the 
present invention: 

(1) Method and Parameter Prediction System for 
Reestablishing a Temporarily Interrupted Dynamic 
Communication Link, U.K. 9613529.8. 



FIELD OF THE INVENTION 

This invention relates generally to 
telecommunications, and in particular, to established 
temporarily- interrupted communication links in an 
dynamic telecommunication system where communication 
parameters such as Doppler frequency offset and 
propagation delay are time- variant . 

BACKGROUND OF THE INVENTION 

In telecommunication systems, links are established 
between nodes of a system for conveying information. 
Efficient use of these links requires enforcement of 
communication parameters defining when nodes at each end 
of a. link anticipate information from the other unit. 
These parameters, in efficient bandwidth systems, comply 
with extremely tight tolerances to maximize information 
exchange. These parameters may include propagation 
delay between nodes of a network and transmit 
frequencies combined with any path introduced frequency 
error such as Doppler frequency shifts. As the 
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communication parameters of the interfacing nodes become 
more predictable and precise, guardbands on 
communication channels such as frequency and timeslots 
can be reduced to accommodate exchange of additional 
5 information . 

In conventional static or lightly dynamic 
telecommunication networks, these communication 
parameters may be refined upon the establishment of a 
communication link and relied upon throughout a 

10 communication session. However, in a dynamic 
communication system where these communication 
parameters are rapidly changing, communicating nodes 
such as subscriber units must continually refine and 
update these communication parameters in order to 

15 maintain the prescribed tolerances of the system. 

Many systems evaluate each transmission employing 
dynamic communication parameters to determine the 
precision of the employed parameters. When the 
parameters exceed a threshold value as defined by system 

20 specifications such as channelized communication 

frequencies or timeslots, the receiving system informs 
the transmitting node of correction values necessary to 
bring the communication parameters into compliance with 
the system requirements. 

25 This approach is suitable when established 

communication links are not susceptible to transient 
effects such as temporary outages due to obstructions or 
fleeting interference. However, in dynamic 
communication systems where communication parameters are 

30 short-lived, even brief outages of communication links 
cause communication parameters to become stale and 
ineffective when utilized following communication link 
outages. A node such as a subscriber unit that employs 
out -dated communication parameters will not be 

35 recognized by a receiving system due to the precise 

tolerances of channelized communication systems. When a 
subscriber unit is not recognized, the communication 
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link will be terminated. Reestablishing a communication 
link is both inefficient and annoying for communication 
link users. 

Furthermore, systems that attempt to predict 
communication parameters based on historical values, may 
be ineffective when a communication link is interrupted 
and one of the communication nodes such as a subscriber 
unit unpredictably alters its direction of travel 
causing an introduction of unanticipated Doppler 
frequency and timing variations. The prediction of 
communication parameters then becomes ineffective when 
the link obstruction ceases. 

Accordingly, there is a significant need for a 
method and apparatus for predicting dynamic 
communication parameters continuously throughout a brief 
outage such that when the outage is terminated, a 
subscriber unit will have revised its dynamic 
communication parameters to comply with the current 
system dynamics. The subscriber unit may then 
immediately employ these revised dynamic communication 
parameters in subsequent transmissions to a receiving 
unit. The receiving unit will recognize and process the 
received communications since the revised communication 
parameters are in compliance with system parameters, 
thus foregoing the need to establish a subsequent 
communication link. 




BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 depicts a highly simplified diagram of a 
5 dynamic communication system of which the present 
invention may form a portion thereof; 

FIG. 2 depicts a block diagram of mobile subscriber 
unit in a dynamic communication system of which the 
present invention may form a portion thereof; 
10 FIG. 3 depicts a block diagram of a communication 

parameter predictor for a subscriber unit of which the 
present invention may form a portion thereof; 

FIG. 4 depicts a block diagram of a vector filter 
for a communication parameter predictor of which the 
15 present invention may form a portion thereof; 

FIG. 5 depicts a block diagram of an acceleration 
compensator in accordance with a preferred embodiment of 
the present invention; and 

FIG. 6 is a flowchart for predicting dynamic 
20 communication parameters in a dynamic communication 

system, in accordance with a preferred embodiment of the 
present invention. 

25 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention has utility in immediately 
resuming communication services upon interruption of a 
communication link between a mobile base station and a 

30 mobile subscriber unit in a communication system wherein 
the mobile base station and the mobile subscriber unit 
are dynamically oriented. In a satellite communication 
system, the dynamic orientation would be due to the 
orbiting nature of the satellite in relationship to the 

35 Earth-bound or near Earth-bound subscriber unit. In the 
preferred embodiment, the mobile subscriber unit is 
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^ highly mobile such as is characteristic of a subscriber 

unit located within an aircraft. 

FIG. 1 depicts a highly simplified diagram of a 
dynamic communication system of which the present 
5 invention may form a portion thereof. A dynamic 

communication system 90 has components or nodes that are 
dynamically oriented such that distances and velocities 
are changing between the nodes. In FIG. 1, dynamic 
co mm u n ication system 90 is shown as an orbiting 
10 satellite -based system where satellites 100 represent 
one of the communication nodes or mobile base stations 
and in the preferred embodiment are space-born. It is 
also practical that dynamic communication system 90 
could have air -born nodes or even that nodes could be 
15 Earth-bound but still in motion in relationship to other 
nodes. Mobile subscriber unit 200 forms the other end 
of communication link 105 in dynamic communication 
system 90. 

The motion or dynamics between nodes creates the 

20 dynamic communication parameters necessary for effective 
and efficient communication in an efficient system. 
Satellite 100 orbits the Earth in satellite orbit 110. 
The term "satellite 11 is intended to include orbiting 
satellites such as low-Earth orbiting (LEO) satellites, 

25 medium-Earth orbiting (MEO) satellites, or highly- 
inclined orbit (HIO) satellites. 

Satellite 100 in satellite orbit 110 orbits the 
Earth at an altitude of around 780 km in a LEO system. 
This low- Earth orbit translates into a velocity with 

30 respect to the Earth of around 25,000 km/hr. Such a 
velocity allows satellite 100 to be in view of a 
terrestrial node such as mobile subscriber unit 200 for 
a duration of only around nine minutes. This velocity 
creates very dynamic communication parameters such as 

35 Doppler frequency shift and propagation delay or timing. 
With such a variation on communication parameters, 
mobile subscriber unit 200 must frequently revise its 
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comr nication parameters to comply with the system 
tinu. 3 as expected by satellite 100. 

In order to establish communication link 105, mobile 
subscriber unit 200 uses default or initial 
5 communication parameters. These parameters may be 
either derived from observations of the broadcast 
channel, pre-programmed into mobile subscriber unit 200 
or may be received by mobile subscriber unit 200 on an 
independent broadcast channel 102. Such a channel is a 

10 broadband channel that is uni- directional from satellite 
100. Broadcast channels may carry other system managing 
information broadcast to all mobile subscriber units 
within view of satellite 100. 

Mobile subscriber unit 200 uses default 

15 communication parameters as an approximation to actual 
communication parameters. This transmission to 
satellite 100 occurs generally on a secondary or 
acquisition channel having a much lower bandwidth and 
broader tolerances for both frequency and time 

20 variations. Satellite 100 evaluates the transmission by 
mobile subscriber unit 200 which employed default 
communication parameters. Satellite 100, in a 
subsequent transmission to mobile subscriber unit 200, 
returns the result of this evaluation to mobile 

25 subscriber unit 200 in the form of feedback 

communication parameters denoting necessary adjustments 
to the default communication pcurameters. In the 
preferred embodiment/ these feedback communication 
parameters are frequency -of -arrival (FOA) raw parameter 

30 225 (FIG. 2) and time-of -arrival (TOA) raw parameter 230 
(FIG. 2). These parameters denote a differential value 
of the expected frequency and timing received at 
satellite 100 and the actual measured values at 
satellite 100. Mobile subscriber unit 200 revises the 

35 default communication parameters according to the 

feedback information. Mobile subscriber unit 200 in 
subsequent transmissions to satellite 100 employs the 
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revised communication parameters. This refining process 
repeats until dynamic communication parameters have been 
sufficiently refined to enable non- interfering 
transmission on a narrow-band channel. Satellite 100 
5 then directs mobile subscriber unit 200 to a traffic 
channel and establishment of communication link 105 is 
complete. This process continues to be iteratively " 
performed as necessary to maintain the communication 
link parameters on the narrow band traffic channel. 

10 FIG. 1 also depicts satellite 100 in motion with 

respect to mobile subscriber unit 200 and the potential 
for an outage of communication link 105. At a first 
point in time, satellite 100 and mobile subscriber unit 
200 are shown to have an active established 

15 communication link 105. 

As mobile subscriber unit 200/ traveling in mobile 
subscriber unit direction of travel 125 , obstructs 
communication link 105' either by banking thus causing a 
aircraft wing obstruction as in the case of an agile 

20 aircraft or passes by an obstruction such as a mountain, 
communication link 105' experiences an outage with 
mobile subscriber unit 200* . For example, in airborne 
applications, a link outage may occur due to control 
surface blockage of the link while the aircraft performs 

25 a maneuver. After the aircraft completes the maneuver, 
the line of sight to the link is reestablished. Unlike 
the ground mobile subscriber cases, the aircraft 
dynamics are not necessarily negligible compared to the 
satellite dynamics. 

30 Also, the aircraft dynamics are time varying during 

the outage time and simple constant parameters are not 
accurate enough for link prediction. For example, as an 
aircraft executes a turn, the aircraft first rolls and 
then the turn gradually ensues. Early in the roll, the 

35 link may become blocked by a control surface, i.e. by 
the wing. If the filter locked- in estimates of the 
dynamics at this point, it would under estimate the 
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dynamics* of the actual turn and not correctly predict 
the link parameters . In order to correctly predict the 
link parameters, the loop tracking filter must account 
for the radial acceleration of the link which comprises 
5 the satellite and the aircraft accelerations. 

When mobile subscriber unit 200' ' emerges from or 
ceases to create an obstruction, both satellite 100'', 
in satellite orbit 110, and mobile subscriber unit 
200", in direction of travel 125 have substantially 

10 changed their orientations with respect to each other. 
If an obstruction obscures communication link 105 for 
even a second in a narrow-band LEO satellite 
communication system, employment by mobile subscriber 
unit 200 9 • of the dynamic communication parameters 

15 utilized at the locations of satellite 100 and mobile 
subscriber unit 200, would be out of the tolerances 
required by communication link 105' 1 . The need for 
revising and predicting dynamic communication parameters 
becomes apparent as mobile subscriber unit 200 operates 

20 in a dynamic environment such as onboard aircraft. 

FIG. 2 depicts a block diagram of mobile subscriber 
unit 200 in a dynamic communication system 90 (FIG. 1) 
of which the present invention may form a portion 
thereof. A receiver 215 receives, through antenna 205, 

25 broadcast channel 102 (FIG. 1) which may explicitly 

contain default communication parameters or be used to 
derive communication parameters. Receiver 215 buffers 
all received transmissions from satellite 100, such as 
broadcast channel information, acquisition channel 

30 information, and information on communication link 105. 

A controller 220 receives the buffered information 
from receiver 215. Controller 220 parses this 
information to extract default or feedback communication 
parameters from the received signal- In a preferred 

35 embodiment, dynamic communication parameters are 
comprised of frequency and timing or propagation 
parameters. Controller 220 sends FOA raw parameter 225 
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and TOA .raw parameter 230 to a communication parameter 
predictor 240. Controller 220 also evaluates the status 
of communication link 105 (FIG. 1) and generates a link 
detect status 235 denoting whether an outage is 
5 underway . 

In the preferred embodiment where mobile subscriber 
unit 200 is located on an aircraft, aircraft avionics 
290 generates acceleration data 291 for use by 
acceleration compensator 292 during a link outage. Link 

10 detect status 235 also enables acceleration data 291 to 
be processed in the vector filter for generation of 
predicted communication parameters during a link outage. 

During the tracking process for maintaining a 
communication link/ raw feedback communication 

15 parameters are employed in transmissions to "pull" 
timing and frequency of mobile subscriber unit 
transmissions into compliance with tight tolerances of 
the traffic channels of communication link 105. The use 
of raw feedback parameters allows subscriber unit 200 to 

20 quickly adjust communication parameters in gross steps 
as compared with finer predicted incremental values 
capable from the vector filter of the present invention. 

Communication parameter predictor 240 receives the 
parameter and status signals and processes them to 

25 generate predicted frequency parameter 245 and predicted 
timing parameter 250. These signals inform frequency 
synthesizer 255 and timer 260 of the requisite 
transmission characteristics required to have 
information transmitted by transmitter 265 of subscriber 

30 unit 200 on communication link 105 received at satellite 
100 during the proper channel {i.e., frequency and 
timeslot in a channelized system) . 

FIG. 3 depicts a block diagram of a communication 
parameter predictor for a mobile subscriber unit of 

35 which the present invention may form a portion thereof. 
Communication parameter predictor 240 improves the 
performance of mobile subscriber unit 200 by providing a 
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smoothing and predicting vector filter 310. Vector 
filter 310 acts upon range measurements (R(k)) 330 and 
velocity measurements (V(k)) 325 and produces smoothed 
estimates of the range and velocity of satellite 100 
5 relative to subscriber unit 200. 

Frequency -of -arrival raw parameter 225 and tiroe-of- 
arrival raw parameter 230 are converted to range and* 
velocity measurements in communication parameter 
predictor 240 for use by vector filter 310. The speed 

10 of light, C/ and a subscriber unit base time defining a 
reference from which subscriber unit transmission and 
reception timing are based/ are used in converting time- 
of -arrival parameters to and from range parameters. 
Smoothed and predicted range estimate 380 and predicted 

15 velocity estimate 375 are converted back to predicted 
frequency parameter 245 and predicted timing parameter 
250 for use by frequency synthesizer 255 (FIG. 2) and 
timer 260 (FIG. 2). 

FIG. 3 also shows the bypass switching of vector 

20 filter 310 as controlled by synchronization complete 
status 238 during the synchronization "pull-in" phase, 
as described above in FIG. 2, for establishing 
communication link 105 (FIG. 1). During 
synchronization/ communication parameter predictor 240 

25 operates with vector filter 310 out of the loop. At 
this point in the operation, the use of a smoothing 
filter is pointless. After "pull-in" is achieved and 
the mode is switched from sync to traffic, vector filter 
310 is switched on provided sufficient raw measurements 

30 have been made to initialize the state vector of the 
vector filter. Initialization of vector filter 310 is 
discussed in FIG. 4. The structure of communication 
parameter predictor 240 converts TOA raw parameter 230 
(FIG. 3) and FOA raw parameter 225 (FIG. 3) into a raw 

35 range parameter 330 and a raw velocity parameter 325 and 
both are input to vector filter 310 as a 2x1 matrix. 
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FIG. 4 depicts a block diagram of a vector filter 
for a communication parameter predictor 240 of which the 
present invention may form a portion thereof. A 
preferred embodiment of vector filter 310 is shown as a 
5 three-state filter. Other numbers of states may be 
implemented for increased accuracy. In vector filter 
310, acceleration is the third state and the dynamically 
random process is now the rate of change of 
acceleration. Other filters that could be considered to 
10 further improve performance would be four or more states 
or and extended Kalman filter. 

The design of a filter requires both a model for the 
physical process and a model for the observation 
process. The physical process model is selected to 
15 incorporate as much information as possible about the 
relationships between range, range-rate, range-rate- 
rate, etc. The more sophisticated the model/ the more 
complex the filter and the better the tracking results. 
The observation model is based upon the measurements 
20 available to the system and the noise process that is 
corrupting these measurements. 

For a vector filter, the process model expressed as 
a first order vector dynamic equation is: 
x(k+l) » A x(k) + w(k) 
25 where w(k) is the physical model noise process. Its 

covariance matrix is defined by: 
Q(k) E{w(k) wT(k)} 
The observation model is given by: 
y(k) = C x(k) + z(k), 
30 where z(k) is the measurement noise process. The 

measurement noise covariance matrix is defined by: 
R(k) = E{z(k) z T (k)J 
The Kalman filter estimator of x(k) given its 
previous value, x(k-l), is given by: 
35 x(k) - A x(k-l) + K(k)[y(k) - C A x(k-l)] 

where the gain matrix is: 



-11- 



BNSOOCIO: <GB 2318948A_I_> 



o 



10 



15 



20 



25 



K(k> = Pl(k) C^lCP^k)^ 1 + R(k)]" 1 
where , 

P 1 (k) « A P(k-l) AT + Q(k-l) 
P(k) - P x (k) - K(k) C(k) P x (k) 

The discrete time difference equation for 
acceleration is found by expressing the acceleration at 
time (k+l)T as the acceleration at time kT plus T times 
the acceleration rate of change. Mathematically the 
radial acceleration, a(k) is: 
a(k+l) = a(k) + Ta' (k) 

where, 

a'{k) is the rate of change of radial acceleration 
at time kT. 

Integrating acceleration to find velocity, v(k): 

^2 



v(k+l) = v(k) + T a(k) + -r- 



a'(k) 



Integrating again to find slant range, r(k) : 



r(k+l) = r(k) + T v(k) + -j- 



•r 3 

a(k) + -g- 



a'(k) 

Defining, 



x(k) 



fr(k)-| 
v(k) 
U(k)J 



x x (k) 
x 2 (k) 



x-^k+l) = x-^k) + T x 2 (k) + -j- 



a'(k) 



x 2 (k+l) - x 2 (k) + T x 3 (k> + y 
x 3 (k+l) = x 3 (k) + T a'(k) 
Expressed in matrix notation: 

*Ti2 



x 3 (k) 



a'(k) 



T3 
6 



x(k+l) = 



1T- 

0 1 T 
LO 0 1 



x<k) 



II 

2 

LT -I 



a*(k) 



which is the first order vector equation. Hence, 
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1T T 

0 1 T 
0 0 1-1 



and 



w(k) = 



T 

2 

l_T J 



a'(k) 



Since independent observations of range and velocity 
are made via the time and frequency tracking loops, 
respectively, the observation equation is: 
5 y(k) = [ J J oJ x <*> + z < k > 

which defines the 2x3 C -matrix for this 
implementation , 

c = r 1 0 °i 

c Lo l oJ 

10 As mentioned above, the structure of communication 

parameter predictor 240 converts TOA raw parameter 230 
(FIG. 3) and FOA raw parameter 225 (FIG. 3) into a raw 
range parameter 330 and a raw velocity parameter 325 and 
are input to vector filter 310 as a 2x1 matrix. A 

15 feedback path is added to the input and multiplied by 
the 3x2 gain matrix K. The result is a 3x1 matrix 
denoted by D. The 3x1 matrix B is added to D to yield 
the output 3x1 matrix which is the smoothed estimate of 
range, velocity and acceleration. The B and D matrices 

20 are storage medium and do not appear directly in the 
vector filter equations. 

Communication parameter predictor 240 and its 
inclusive vector filter 310 must operate through periods 
of time when no data is available for input but an 

25 updated estimate is still needed for the tracking loop. 
This occurs routinely when coasting in-between 
maintenance bursts containing feedback communication 
parameters or sporadically when communication link 105 
blockages occur. The availability of a measurement can 

30 be determined by link detect status 235. When 

available, vector filter 310 operates closed loop and 
all covariance matrices are updated. When no data is 
available, the feedback tap from the loop gain matrix is 
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opened/ the gain matrix and covariance matrices are not 
updated, and the estimate is formed from the geometric 
matrix A, the previous state estimate matrix, and 
acceleration compensation data from aircraft avionics 
5 such as an aircraft inertial navigation system. 

Since the state variables are updated at the frame 
rate and the gain and covariance matrices are updated at 
the maintenance schedule rate which is every fourth 
frame in the preferred embodiment, two A matrices must 
10 be maintained. The matrix for the state variable 
calculations is designated as Al and is given by: 
1 Ti Ti 2 /2 



Al = 



0 1 Ti 
0 0 1 



The matrix for the gain and covariance matrix 
calculations is computed for a four frame interval. It 
15 is designated as matrix A4 and is given by: 

1 T4 T4 2 /2 



A4 = 



0 1 T4 
0 0 1 



Cain Matrix and E rror Covariance Matrix Computation 
The procedure for computing the error covariance 

20 matrices, P and PI, of an error covariance matrix 

computation block 410 is discussed below. The first 
time the loop is executed, an initial value for the P 
matrix, P(0), must be available. The initialization 
values used are discussed in a later section. From 

25 this, Pl(l) is computed. The gain matrix, K(l), of a 
gain matrix computation block 420 is then computed and 
used in the filter loop of vector filter 310 to update 
the state variable estimates. The gain matrix is also 
used to update the P matrix, P(l), which is stored until 

30 the next time the loop is executed. 

Upon subsequent executions of the loop, the previous 
value of the P(k) is used in the Pi(k) and K(k) 

computation rather than the initial value matrix. 
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However., if it is determined by subscriber unit 200 that 
a link outage has occurred then the P matrix will be set 
back to P(0). This allows reinitialization of the error 
variances back to their initial values and increases the 
loop gain to facilitate rapid reacquisition of the 
uplink from subscriber unit 200 to satellite 100. 

Mobile subscriber unit 200 should also reset the* 
state vector to initial values after a link outage. 
Initial values are determined by feedback communication 
parameters after communication link 105 is re- 
established. 

A simplification of the P matrix computation is 
shown implemented in error covariance matrix computation 
block 410 and is defined by substituting the values for 
the C matrix and factoring Pi from the equation: 

P(M - Pl(k) - K(k) C(k) Pi(k) 

- [I - K(k) C(k) ] Pi(k) 
-l-Kn(k) Ki2(k) 0- 

K 2 l(k) 1-K 2 2(k) 0 p 1(Jc) 
L K 3 l(k) K 32 (k) lJ 



As the P matrix is iterated over and over again, its 
entries will stabilize to their steady state values. 
The entries in the P matrix change from their initial 
error variances of the estimates to their final values. 

A gain matrix computation block 420 is shown in FIG. 
4. The K matrix is computed from the PI matrix and the 
other constants C and R. The dimensionality of each 
matrix used in the computation is also shown. Note that 
even though a 3 -state filter is implemented in the 
preferred embodiment, the inverse function acts on only 
a 2x2 matrix. As vector filter 310 is iterated the 
filter gains approach their steady state values. 

Measurement Noise Covarian^ Matrix 

The measurement noise, z(k), is defined as: 
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f"zi(k)~j _ r Range Measurement Noise 1 

2 ^ ^ =5 LZ2(k)J ~ L Velocity Measurement NoiseJ 

The range and velocity measurements are derived 
5 directly from the TOA raw parameter 230 (FIG. 3) and FOA 
raw parameter 225 (FIG. 3) respectively. Since the . 
variances of TOA raw parameter 230 and FOA raw parameter 
225 are known, the range and velocity variances can be 
determined. The measurement processes are mean zero 
10 (i.e., no bias error), and hence the variance = E{ 

zi 2 (k) } which is what is needed to define terms in the 

R matrix. 

In the preferred embodiment, FOA raw parameter 225 
standard deviation is: 
15 Sf - 31.97 Hz at 4.1 dB Eb/No 

Hence , 

c 3 x 10 8 _ Q 

s v = sf —r- = 31.97 T =5.9 

v 1 fc 1625 x 10 6 

m/s 

20 and, 

s v 2 * 34.84 m 2 /s 2 

Likewise, TOA raw parameter 230 standard deviation 

is: 

25 st - 0.0028 ms at 4.1 dB Eb/No 

Hence, 

s r = c st = 3 x 10 8 x 0.0028 x 10" 3 = 840.0 m 
and, 

30 s r 2 = 705,600 m 2 

As for the off diagonal terms, these terms are the 
cross -correlations between the range and velocity 
measurements. These measurement errors are considered 
35 uncorrelated for the following reasons: 
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(i) • Under static conditions, TOA raw parameter 230 
is a more robust measurement than FOA raw parameter 225. 
It has a very small standard deviation over it's 
operating signal -to -noise ratio where FOA raw parameter 
225 has a much broader range; 

(ii) Under real orbital conditions, TOA raw 
parameter 230 and FOA raw parameter 225 processes arfe 
not stressed at the same time. TOA raw parameter 230 
has large error values at low elevation angles where 
radial velocity is maximum and acceleration is minimum, 
and FOA raw parameter 225 has large error values at high 
elevation angles where radial acceleration is maximum 
and velocity is minimum; and 

(iii) Under link blockage or fading conditions, no 
FOA raw parameter 225 and TOA raw parameter 230 updates 
are made to vector filter 310 (or the tracking loop for 
that matter) so correlated error under this condition is 
not a consideration. 

The R matrix for this filter implementation is: 



Process Noise Hovarianne Matriy 

Filter development relies on the modeling of radial 
acceleration. The true behavior of the acceleration is 
that it is near zero when satellite 100 is viewed at low 
elevation angle and is some maximum negative value at 
zenith, satellite 100 has slowed down relative to mobile 
subscriber unit 200. The magnitude of the acceleration 
then reduces back to zero as the elevation angle 
decreases back to zero. Non-overhead passes have a 
similar profile without the maximum acceleration 
magnitude being reached. The acceleration from one 
transmission from mobile subscriber unit 200 to the next 



R(k) = E{z(k) z T (k)} 




E{zi 2 (k)} E{zi(k)z2(k)} 
z 2 (k)zi(k)} E{z 2 2 (k)} 
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will have significant correlation as the acceleration is 
a slow varying process. However, vehicular motion of 
mobile subscriber unit 200 may introduce other 
extraneous accelerations which may reduce the 
5 correlation from transmission to transmission. 

In the preferred embodiment, the acceleration rate- 
of -change is characterized with values of ±0.5 m/s 3 for 
the overhead pass and ±0.3 m/s 3 for the side pass 
utilized. Maximum acceleration occurs when satellite 

10 100 is directly overhead of subscriber unit 200, In the 
preferred embodiment/ the acceleration rate of change is 
maximum about a minute before and after satellite 100 is 
directly overhead subscriber unit 200 and has a 
magnitude of about 0 . 5 m/s 3 . This corresponds to about 

15 a 70° elevation angle. 

A histogram, Pa'(a')> of the acceleration rate-of- 
change may be constructed to evaluate the variance of 
the process. 



E{a f (k) 



a'(k)} « J a'2 p a '(a') da 1 - 0.053 



20 m 2 /s 3 

Knowing the acceleration rate -of -change variance in 
the preferred embodiment, the covariance matrix can now 
be evaluated. 

Q(k) = E{w(k) wT(k)} 



T T 3 /6 ] 



25 = E{ T 2 /2 a'(k) [ T 3 /6 T 2 /2 T] a'(k) 



[T 6 /36 T 5 /12 T 4 /6 
T 5 /12 T 4 /4 T 3 /2 
T*/6 T 3 /2 T 2 . 



E{a'(k) a'(k)} 



setting T=0.36 sec (the feedback parameter loop 
30 update rate of the preferred embodiment) and using the 
result for the expected value of a' 2 , the Q matrix has a 
value of: 
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25 



T3.5: 
2.71 
Ll .41 



.51e-6 2.67e-5 
Q(k) « |2.70e-5 2.23e-4 
l8e-4 1.24e-3 




Filter Initialization 
5 Initial values are needed for the filter state 

variables (r(0), v(0) and a(0)) and the 3x3 prediction 
mean square error matrix P(0). 

The state variables r(0) and v(0) are initialized by 
the range and velocity estimates determined during when 

10 subscriber unit 200 is performing a synchronization mode 
while it is "pulling in" to meet the criterion of 
communication parameters for communication link 105. 
a(0) is initialized by differencing two velocity 
measurements; [ v(l) - v(0) ] / T . Note that to 

15 initialize a(k), two measurements are needed* Since 

communication link 105 is now operating in Traffic mode 
and hence, in the preferred embodiment, these 
measurements are about 0,36 sec apart. In subscriber 
unit 200, tracking loops will run without vector filter 

20 310 until these measurements are available to the 
filter. 

The mean square error matrix P(k) should be 
initialized by the variance of the initial state 
measurements. The P(0) matrix is given by:. 



fE{r*r 
P(0) = E{vr 
LE{a*r 



} E{r-v] 
} E{v»v} 
] E{a»v} 




where f 

E{r»r} = s r 2 

E{r»v} = E{v»r} = 0 

30 E{r»a} » E{a»r} = 0 

E(v»v} = s v ^ 

vi -vn Sv^ 

E{va} = E{a.v} = E{ ^ U vU = 

and, 
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vi-vo vi -vp 
E{a»a} = E{ — £ — — £ — 



_ 2Sy 2 



10 



Updating the matrix: 



P(0) = 



s r < 
0 

0 



Sv 2 
Sy^ 



Sy^ 
T 

2Sy 2 



T T 2 

Substituting the values given in the last section 
for s r 2 and s v 2 and using T=0.36 sec since the initial 
values are calculated over a four frame interval, the 
initial value, in the preferred embodiment, for the 
error covariance matrix is: 



t 705600 0 
0 34.84 96 
0 96.78 537 



° 1 
• 78 

7.65J 



FIG. 5 depicts a block diagram of an acceleration 
compensator in accordance with a preferred embodiment of 

15 the present invention. When a link interruption occurs, 
the loop tracking filter combines the geometric matrix 
A, the previous state estimate matrix, and acceleration 
compensation data from aircraft avionics. In order to 
correctly predict the link parameters, the loop tracking 

20 filter must account for the radial acceleration of the 
link which comprises the satellite and the aircraft 
accelerations. The tracking loop filter calculates the 
radial range, velocity, and acceleration based on the 
observations made of Doppler offset and propagation 

25 delay. These radial parameters have components due to 
the aircraft dynamics as well as the space vehicle 
components. The radial acceleration can be expressed 
as : 

30 a=a sv • a A 

where , 
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a is the scalar value of the composite radial 
acceleration; 
— > 

a sv is the satellite acceleration vector; 
5 -> 

a A is the aircraft acceleration vector; and 
• denotes the vector dot product. 

When the link outage occurs and no observations of 

10 the satellite are available, the last estimate of the 
satellite component of radial acceleration should be 
held constant and the aircraft portion should be 
continually updated. However, in order to do this, the 
two components of the radial acceleration must be 

15 isolated . 

In airborne applications which use a gimbaled dish 
antenna for geosynchronous satellite communications, a 
bore sight accelerometer can be mounted to the antenna. 
The bores ight accelerometer gives a continuous reading 

20 of the aircraft acceleration in the direction of the 
satellite with a significant satellite component. 
However, in many airborne applications, gimbaled 
antennas are replaced by a patch or fin type antennas 
and bore sight accelerometers are not desirable since 

25 spatial tracking is lost during a link blockage and the 
output of the accelerometer is inaccurate during an 
outage . 

In the preferred embodiment, an aircraft 
accelerometer vector is derived from the aircraft 

30 inert ial navigation system combined with the link 
observations. The composite scalar acceleration is 
calculated by the Kalman filter based on observations of 
the link Doppler and propagation delay. The aircraft 
acceleration vector may be calculated from the velocity, 

35 roll, pitch, yaw, and heading information available from 
the aircraft inert ial navigation system. The satellite 
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acceleration vector is then calculated by dividing the 
scalar observation by the aircraft acceleration vector. 

— > -» 

a sv « a/a A 

5 

All three components are know known and when an link 
outage occurs, the satellite term is held constant at 
the last estimate and the aircraft term is continuously 
updated. The computation needed to perform this 

10 calculation is shown in FIG. 5. 

FIG. 6 is a flowchart for predicting dynamic 
communication parameters in a dynamic communication 
system, in accordance with a preferred embodiment of the 
present invention. In the preferred embodiment, mobile 

15 subscriber unit 200 performs link outage protection 
procedure as shown. 

In a task 504, mobile subscriber unit receives 
acceleration data from an aircraft inertial navigation 
system. The mobile subscriber unit, in the preferred 

20 embodiment, is located on an aircraft with the aircraft 
capable of generating navigational data from which 
acceleration may be calculated. In a task 505, mobile 
subscriber unit 200 initiates the establishment of 
communications with satellite 100 by initiating 

25 communications as directed by a user of mobile 

subscriber unit 200 either in response to notification 
that an incoming call is pending or as initiated by the 
user. 

A task 510 then receives feedback communication 
30 parameters from satellite 100 either from an independent 
broadcast channel or in a maintenance portion of 
communication link 105. These feedback communication 
parameters , in the preferred embodiment, comprise 
frequency and timing variances associated with 
35 transmission of data from mobile subscriber unit 200 to 
satellite 100 as evaluated by satellite 100. 
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^ A query task 515 then evaluates if synchronization 

is completed. Synchronization is complete when 
communication parameters are sufficiently accurate that 
satellite 100 permits mobile subscriber unit 200 to 
5 communicate on a tighter- tolerance traffic channel. 
Also, synchronization is complete, in the preferred 
embodiment/ when the matrices of vector filter 310 are 
initialized. 

A task 550 proceeds when the learning or 

10 synchronization of query task 515 is not complete. In 
task 550 , the feedback or raw parameters of 
communication link 105 are employed for communicating 
with satellite 100. These feedback parameters provide 
larger "pull- in" increments while coarse adjustment of 

15 communication parameters is underway. 

When synchronization is completed, a task 520 
utilizes the communication parameters emerging from 
vector filter 310. These parameters more incrementally 
continuous and provide a smoothing effect to 

20 communication parameters employed on communication link 
105. This smoothing provides a collateral benefit to 
communication link 105 in that vector filter 310 tracks 
variations in communication parameters and anticipates 
future changes to these parameters. This anticipation 

25 and employment of these anticipated communication 
parameters reduces the frequency of transmission of 
feedback parameters on communication link 105. 

A task 521 calculates the acceleration vector for 
the satellite. When a link outage occurs and no 

30 observation of the satellite are available, the last 
estimate of the satellite component of radial 
acceleration is saved and held constant while the 
aircraft portion is continually updated during the 
outage . 

35 A query task 530 evaluates communication link 105 to 

determine if the link remains intact or if an outage has 
occurred due to link blockage. If link interruption has 
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i) 
A 

not occurred/ then processing returns to task 510 to 
listen for any subsequent transmissions of feedback 
parameters to be incorporated into future predictions of 
communication parameters . 
5 When an link outage is detected/ a task 531 

retrieves current acceleration data for computation with 
predicted values as discussed in FIG. 5. A task 535' 
predicts future communication parameters for use in 
communication link 105. These parameters are generated 
10 when link detect status 235 disables the input path and 
vector filter 310 cycles producing predictions of 
communication parameters. 

A task 540 then employs these predicted 
communication parameters by attempting to transmit to 
15 satellite 100 using these predicted communication 

parameters. These predicted parameters roughly track 
the actual communication parameters necessary to 
reestablish communications . 

A query task 545 then reevaluates communication link 
20 105 to determine if the link blockage has ceased. If 
the link has not been successfully reestablished, then 
processing returns to task 535 for a subsequent 
prediction of the communication parameters presently 
necessary to attempt a successive re -establishment of 
25 communication link 105. This iterative prediction of 
subsequent parameters may continue indefinitely until 
the link is recovered or may be terminated after a 
predetermined number of attempts. Because the 
prediction of subsequent communication parameters has 
30 inherent error which Depends on the number of stages 
utilized in vector filter 310, the prediction of 
communication parameters will drift from actual 
necessary values for the re -establishment of 
communication link 105. A predetermined timeout or an 
35 automatic re- initiation of communications as described 
in task 505 may be employed. 
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If -re -establishment of communication link 105 is 
detected in query task 545, then processing passes back 
to task 510 for reception of feedback communication 
parameters and reinitialization of vector filter 310. 
5 The present invention has utility in maintaining 

synchronization with a satellite and a mobile subscriber 
unit during link blockages such as those due to control 
surface blockage of the link while the aircraft performs 
maneuver . 

10 The present invention has been described above with 

reference to a preferred embodiment. However, those 
skilled in the art will recognize that changes and 
modifications may be made in this preferred embodiment 
without departing from the scope of the present 

15 invention . 
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CLAIMS 

Waat is claimed is: 

1. In a wireless communication system, a method of 
5 restoring a communication link between a mobile base 
station and a mobile subscriber unit following an 
interruption of said communication link comprising the 
steps of: 

establishing said communication link having dynamic 
10 communication parameters; 

revising said dynamic communication parameters as 
orientation changes between said mobile base station and 
said mobile subscriber unit; 

upon an interruption of said communication link, 
15 providing acceleration data characteristic of said 
mobile subscriber unit; 

said mobile subscriber unit predicting current 
dynamic communication parameters in part from said 
acceleration data; and 
20 upon termination of said interruption of said 

communication link, said mobile subscriber unit 
employing said current dynamic communication parameters 
for restoring said communication link. 

25 2. The method as recited in claim 1 further 

comprising the step of generating said acceleration data 
representative of motion of said mobile subscriber unit. 
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3 . The method as recited in claim 1 wherein said 
revising step comprises the steps of : 

transmitting on said communication link using said 
dynamic communication parameters; 

evaluating at said mobile base station an accuracy 
of said dynamic communication parameters employed by* 
said mobile subscriber unit; 

when said accuracy of said dynamic communication 
parameters exceeds a threshold/ transmitting feedback 
communication parameters denoting adjustments to said 
dynamic communication parameters to said mobile 
subscriber unit; 

receiving said feedback communication parameters; 

and 

employing said feedback communication parameters to 
generate said dynamic communication parameters for use 
in subsequent transmissions by said mobile subscriber 
unit. 

4 . The method as recited in claim 3 wherein said 
employing step further comprises the steps of: 

when said mobile subscriber unit receives feedback 
communication parameters, incorporating said feedback 
communication parameters into said dynamic communication 
parameters ; and 

else generating dynamic communication parameters 
from previous ones of said dynamic communication 
parameters . 

5. The method as recited in claim 4 further 
comprising the step of employing a Kalman-type filter in 
said incorporating and generating steps. 

6. The method as recited in claim 1 wherein said 
predicting step further comprises the steps of: 
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filtering previous ones of said dynamic 
communication parameters; 

incorporating said acceleration data; and 
predicting current dynamic communication parameters . 

7. The method as recited in claim 6 further 
comprising the step of attempting to reestablish said 
communication link using said current dynamic 
communication parameters as said dynamic communication 
parameters - 

8. The method as recited in claim 6 further 
comprising the step of employing a Kalman-type filter 
for performing said filtering steps. 

9. The method as recited in claim 1 further 
comprising the steps of employing Doppler frequency 
offset and propagation delay as said dynamic 
communication parameters . 

10. The method as recited in claim 1 wherein said 
mobile base station is an orbiting satellite and said 
mobile subscriber unit is aboard an aircraft. 
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